The Wiedemann-Franz law in the SU(N) Wolff model 
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We study the electrical and thermal transport through the SU(N) Wolff model with the use of 
bosonization. The Wilson ratio reaches unity as N grows to infinity. The electric conductance is 
dominated by the charge channel, and decreases monotonically with increasing interaction. The 
thermal conductivity enhances in the presence of local Hubbard U . The Wiedemann-Franz law is 
violated, the Lorentz number depends strongly on the interaction parameter, which can be regarded 
as a manifestation of spin-charge separation. 
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The Wiedemann-Franz law is one of the basic prop- 
erties of a Fermi liquid^. Simply stated, it reflect the 
fact, that the ability of quasiparticles to carry charge is 
the same as to transport heat. A possible breakdown of 
this relation is interpreted in terms of spin-charge sep- 
aration. Spinons can carry heat in much the same way 
as charges do, and contribute to heat transport. On the 
other hand, they fail in electric transport, while charge 
excitations contribute to both. Another explanation re- 
quires inelastic scattering. 



This phenomenon is demonstrated in the SU(N) Wolff 
modelSiS*^. It consists of N species of electrons inter- 
acting with each other only at a single site. The model 
is studied with Abelian bosonization^ for arbitrary N. 
We consider fermions with SU(N) spin index or alterna- 
tively we take the ensemble of spin and orbital degrees 
of freedom into account by the N indejsLS,. These ad- 
ditional degrees of freedom can be realized through or- 
bital degeneracy, for example, as in Mn oxides^. As a 
result, the additional degrees of freedom can be called 
flavour or color index. The Wolff model is one of the sim- 
plest impurity models, where electron correlation is still 
presenbSi^iiS. Also it is the basis of studying the effect 
of Coulomb interaction on resonant tunneling through a 
single quantum leveUiii^. 



The Wilson ratio approaches unity as N grows, it 
crosses over to non-interacting (mean-field) behaviour 
as N grows to infinity, a phenomenon characteristic to 
SU(N) modelsi^iii. The electric conduction involves only 
the charge sector, and increasing interaction increases the 
resistivity. Heat is transported by both spin and charge 
excitation, and heat conduction is favoured by any finite 
value of U, regardless to its sign. The respective conduc- 
tivities depend strongly on the local interaction, leading 
to the violation of the Wiedemann-Franz law even in the 
large N limit. This suggests that spin-charge separation 
is the origin of this breakdown. 



The Hamiltonian describing N different species of elec- 



trons interacting only at the origin is given by: 
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and only the radial motion of the particles is accounted 
for by chiral (right moving) fermion fields^^. The model 
can be bosonized via&i& 



j,i\/47r'I>m(a;) 



and after introducing charge and spin fields aa^ 
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rt = 1. . . — 1, the Hamiltonian separates into different 
sectors: the spin sector is described by N— 1 identical de- 
coupled bosonic modes, and the charge sector transforms 
into a similar massless bosonic mode as 
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These can readily be diagonalized by introducing pure 
bosonic representation of the <i>(a;) fields as 
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where a is the ultraviolet cutoff. Usually x S 
[— L/2,L/2], but L ^ oo in actual calculations. The 
charge Hamiltonian is rewritten as 

Hc = Y.\ vqb+b, + ^Ed^ib, ~ b+) + 
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which describes spinless bosons scattered by U, and a 
source term E. The latter can be transformed out by a 
linear shift of the bosonic field, while the former can be 
considered exactly via Dyson equation. Similar equa- 
tions describe the spin sector. Within the realm of 
bosonization, U is restricted to {—Uo/{N — 1), Uq) with 
Uq — ttv/uq, Uq is the average density per spin in the 
homogeneous case^. 

As a first step to understand the response of our model, 
it is instructive to investigate the Wilson ratio. The Wil- 
son ratio characterizes to what extent the impurity in- 
teraction influences the conduction electron properties. 
Using Ref. 5, it is calculated for the SU(N) Wolff model 
as 
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It interpolates smoothly between zero in the attractive 
case to N/{N — 1) in the repulsive case. Interestingly, 
the latter value was calculated in the N-fold degener- 
ate Coqblin-Schrieffer model— as well. This might sug- 
gest that the U ^ Uq limit of the Wolff model found 
by bosonization may describe the U ^ oo limit of the 
original model. The shortcoming of bosonization is to 
underestimate the strength of interaction as it does for 
simple impurity scattering^. 

As N grows, R approaches unity (i.e. the non- 
interacting limit). Indeed, the ground state of the N 
oo limit of impurity models is the mean- field solution^^, 
because fluctuations around the saddle-point are sup- 
pressed by 1/iV, and vanish as N reaches infinity. As 
a result, both quantities determining the Wilson ratio 
are proportional to the density of states at the Fermi 
energySi, hence their ratio is 1. The general behaviour of 
the Wilson ratio for various A^'s is shown in Fig. ^ The 
SU(2) case was also investigated in Ref. llQ. 

The calculation of the conductivity can be carried out 
in several alternative ways. The most common one relies 
on the physical assumption that the whole voltage drop 
occurs at the impurity site. The total number of electrons 
to the right and left of the impurity are A'^^ and N^, 
respectively. Then the current satisfying the continuity 
equation reads asSii^ 

/ - ^dt {Nb. - Nl) = [H, Nr - Nl] . (10) 

From this, the linear conductivity can be calculated by 
the Kubo formula. 
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FIG. 1; (Color online) The Wilson ratio of the SU(N) Wolff 
model is plotted for A = 2, 3, 4, 8 and oo (red) from top to 
bottom. 



On the other hand, one might think, that the conduc- 
tivity at the impurity site is given by the autocorrelator 
of the current flowing through it. It turns out that these 
definitions of the current yield identical result, and we are 
in position to calculate the conductivity. From the con- 
tinuity equation {dtn + dxj = 0), the local conductance 
at t he imp urity site is obtained, if we proceed following 
Ref. Illligj The local density is determined by 



n{x) = \ —dx^{x), 



(11) 



from which the local current operator is calculated as 



j[x) = -e\j —vdx^{x) 
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Only charge excitations participate in electric trans- 
port, spinous only transport heat, not charge, as will be 
demonstrated below. The expectation value of dx^{x) is 
finite as seen in Ref. 0, which is a natural consequence of 
the same chirality of all A^ channels. The current-current 
correlation function can be evaluated from the bosonic 
representation of the <^(x) field through the equation of 
motion method. After solving the closed set of equations 
using Eq. the frequency dependent ac conductance 
is evaluated as 
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where Xo('^) is the local charge correlator per spin a,t U — 
0. Identical result was obtained for the charge correlator 
using the generating functionai^. It is shown in Fig. |21 
for A^ = 4. Similar curves describe the N ^ 4 behaviour 
as well. Its dc part yields to 
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which decreases smoothly for repulsive interaction, but 
increases rapidly in the attractive sector. Throughout 
the calculations we set ?i = 1, which explains the e^/27r 
factor in front of the right hand side of Eq. ^1 By restor- 
ing the original units, this gets replaced by the universal 
channel conductance unit, /h. Qualitatively similar be- 
haviour was identified in the conductance of a one dimen- 
sional interacting electron gas (Luttinger liquid) through 
a weak link^. Repulsive interaction caused perfect re- 
flection, attractive one resulted in perfect transmission. 
Although the details are different, our results exhibit the 
same phenomenon as seen from Eq. 1141 and from the in- 
set of Fig. [5] In the N— oo limit, the dc conductance 
is zero for any finite U . The charge sector involved in 
electric transport is governed by Eq. © , where the scat- 
tering term (cx U) grows with N, hindering any transport 
in the mean-field limit. 

The ac conductance of an Anderson impurity was in- 
vestigated in Ref. 21 in a different regime of the model, 
and hence with different conclusions from ours. 
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FIG. 2: (Color online) The frequency dependent conductance 
is shown for = 4 for U /Uo = -0.33, -0.3, -0.2 below 
the dashed line indicating U = Q with increasing spectral 
weight. The other curves correspond to U /Uo = 0.2, 0.5 and 
0.9 with decreasing peak close to = W , where VF is a sharp 
momentum space cut-off. The inset shows the dc conductance 
as a function of the Hubbard interaction. 

The heat transport can be studied similarly to the elec- 
tric one. The heat-current operator is obtained from the 
en ergy conservation, following the steps outlined in Ref. 
I22l23 . and reads as 
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each other, hence their contributions can be determined 
independently. The thermal conductivity can be cal- 
culated from energy current-current correlation function 
given by 
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and the retarded response function is obtained after an- 
alytic continuation as 

K Im(n(ia'„ — > w -I- i5)) 
T ^ T\} ■ 

Its calculation amounts to determine the 
(r^9^($(0,r))2(9^$(0,0))2) correlator. Possible 
pairings are of {Trdx^{0, t)92.$(0, 0)) type in all sectors, 
which equals to the the charge and spin correlation 
functions, obtained in Ref. 5. The Hamiltonians describe 
scattering of free bosons on a localized " impurity" , hence 
only self energy corrections can be taken into account 
(vertex corrections are absent). The calculation can be 
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FIG. 3: (Color online) The Lorentz number in the SU(N) 
Wolff model is plotted for N = 2 (solid line), 3 (dashed line), 
4 (dashed-dotted line), 8 (dotted line) on a semilogarithmic 
scale. 

carried out in a straightforward manner, and by taking 
all the channels into account, we find 
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where Gc{x) = v^TTXoix)/il + [N - l)Uxoix)) and 
Gs{x) = v'^TTXo{x)/il - Uxo{x)). As T ^ 0, this can 
readily be evaluated to yield 

TT 



The Hamiltonians in the different sectors commute with 
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In the absence of U, it gives the pure result k = {tt/6)NT 
(or Tr^k^NT/3h upon restoring original units). In the 
presence of finite U, it increases regardless to the sign of 
the interaction. In other words, local Coulomb repulsion 
or attraction facilitates heat transport. From these, the 
Lorentz number reads as 
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where Lq = 7r^fc^/3e^ is the non- interacting result. For 
any finite U and N, the Wiedemann-Franz law is broken, 
except a single value of attractive U for any given N. 
This breakdown is a natural consequence of spin-charge 
separation, as seen in Eqs. ^ and ©. Had we chosen 
iV = 1, the Wiedemann-Franz law would hold. As U in- 
creases in the repulsive regime, it enhances strongly. In 
the attractive case, a minimum value is reached for de- 
creasing U before the divergence at —Uo/iN—1). For fi- 
nite T, the Lorentz number decreases monotonically. Vi- 
olation of the Wiedemann-Franz law was reported in one- 
dimensional electron gas (Luttinger liquid) with repul- 
sive interaction. For any finite interaction, L increased 



from the Fermi liquid value, similarly to our case. From 
these similarities in transport and in X-ray response^, we 
conclude that the physical quantities evaluated within 
our model and its lattice version (Hubbard model) ex- 
hibit the same type of behaviour as a function of inter- 
action. Other failures of the Wiedemann-Franz law have 
been seen in mean-field theories as wellSSiSi. It is plotted 
in Fig. |3|for various N's. Although the long time asymp- 
totics of the Green's function is of Fermi liquid character, 
still spin-charge separation occurs on microscopic level, 
as seen in Eqs. (O and ©, causing the breakdown of the 
Wiedemann-Franz law. Similar phenomenon is expected 
to occur in the single- and multichannel realizations of 
the Kondo effect^. 

In conclusion, we have studied the electric and 
heat transport properties of the SU(N) Wolff model. 
As N grows to infinity, the systems crosses over to 
non-interacting behaviour. At the same time, the 
Wiedemann-Franz law remains broken due to spin-charge 
separation. 
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